Fractal to Nonfractal Phase Transition in the Dielectric Breakdown Model 
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A fast method is presented for simulating the dielectric-breakdown model using iterated conformal 
mappings. Numerical results for the dimension and for corrections to scaling are in good agreement 
with the recent RG prediction of an upper critical rjc — 4, at which a transition occurs between 
branching fractal clusters and one-dimensional nonfractal clusters. 



Diffusion- limited aggregation j|](DLA) is one of the 
most interesting and difficult fractal growth models, de- 
spite its simple rules. It models a wide range of processes 
in nature that give rise to branching structures, including 
viscous fingering ||], electrodeposition and dendritic 
growth [Q. Yet, despite the importance of the model 
and various theoretical approaches ||], a full, controlled 
understanding is lacking. 

An extension of DLA is the dielectric-breakdown model 
[^(DBM), also known as the rj model. This model pro- 
vides a continuously varying fractal dimension d, rang- 
ing from 2 (?7 = 0, the Eden model (7|J^), to roughly 
1.7 (77 = 1, DLA), to 1. Recently |||, it has been ar- 
gued that d = 1 for rj > rjc = 4:, leading to a controlled 
renormalization group expansion for the DBM. Further, 
while continuum DLA has certain integrability proper- 
ties, most other theoretical approaches do not rely on 
these and should apply to the DBM as well. Thus, nu- 
merical studies of the DBM provide a strong check on any 
theoretical approach to fractal growth. Given the intense 
recent activity in the field, stimulated by the conformal 
mapping approach such a check is crucial. 

Several numerical studies of the DBM have been per- 
formed |]l^-^16|] , and in one such study |l^ an r/c equal to 
4 was suggested. However, the difficulty of numerically 
solving Laplace's equation has restricted the cluster sizes 
to roughly 5000 particles for 77 = 2 and to roughly 1000 
particles for larger rj. We show below that the numeri- 
cally measured dimension of clusters with 77 > 4 decreases 
with increasing cluster size, and that 1000 particles is not 
sufficient to extract correct information. 

In this work, we present a fast algorithm for the DBM 
using iterated conformal mappings, leading to an order 
of magnitude increase in cluster size. Unlike numerical 
solution of a discrete Laplace equation, this technique is 
free of grid effects and boundary effects. The numerical 
results are in good agreement with upper critical rj found 
in the RG |^], and for rj > rjc, we find that the finite size 
corrections to the numerically measured dimension are 
in good agreement with the corrections predicted by the 
RG, and extrapolate well to D — 1. 

Conformal Mapping Method — The conformal map- 
ping technique for simulating DLA proceeds as follows: 
define a function F^'^\z) that maps the unit circle in the 
complex plane onto the boundary of the growing clus- 
ter after n growth steps. All singularities of F must lie 



within the unit circle, which implies that F^"'\z) has a 
Laurent expansion 

F(")(z) = f[''\ + ^(j") + + ... (1) 

Take F(o)(z) = z. Then, to obtain from F^"', 

pick a random angle 0„+i, and construct an elementary 
mapping which produces a bump of linear size 

\n+i at angle 6n+i ■ Then, compose mappings to define 
F("+i)(z) = (z)), which adds a bump to the 

cluster. To obtain bumps of fixed size in the physical 
plane, we chose 
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«+i - |^(„)'(e»e„+i)|2' (2) 

where Aq is some constant. Various elementary functions 
are possible; we use the function / proposed in Ref. [p^ , 
which includes a parameter a describing the bump shape. 

For the dielectric breakdown model with given 77, we 
need to chose the angle 9n+i with probability propor- 
tional to |F(")'(e*''"+i)|i"''. To do this, the program uses 
a Monte Carlo technique. It starts with an angle 
which is initially taken equal to On- Then, it picks a 
random "trial angle" 9. Define 



P 



|F(n)'(ei9old)| 



(3) 



If P > 1, it replaces 6'oid with 9, "accepting" angle 9. If 
P < 1, it replaces 6'oid with 9 with probability P. The 
program proceeds through a certain number T of these 
trials, and then takes 9n+i = 9o\d to add the n + 1-st 
particle to the cluster. If T is large enough, this yields 
the desired distribution. 

The time for the algorithm to produce a cluster of size 
N scales as N'^T. We have used various means to deter- 
mine the minimum value of T needed. In Fig. 1, we show 
a plot of Pi against number of growth steps for single real- 
izations with ry = 4, iV = 10000 and T = 50, 300, 400, 500. 
While the curve for P = 50 is clearly different, the other 
curves are all close, suggesting that T — 300 suffices. 
Determinations of the dimension using T — 300 are also 
indistinguishable from those using T = 500; however, 
to be careful, in all quantitative results given below, we 
have chosen T — 500, except for rj — 2, where P = 200 
was used. On an 850-Mhz Pentium HI, a 10000 parti- 
cle cluster with P — 500 can be simulated in roughly 33 
hours. 
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To provide an analytical estimate of the minimum T 
needed we consider the worst case, a one-dimensional 
cluster, obtained by a suitably regularized version of the 
mapping F{z) = Fiz + Fx/ z. For 77 > 2, the proba- 
bility distribution \F^'^'> (e*^)|^~'' diverges near 9 = 0,ti. 
The appropriate regularization of the cluster at a physi- 
cal length of order a/A cuts off this divergence at an angle 

that scales as \J VX/Fi cx \f\Jn. In order to have suffi- 
cient trials that the algorithm is able to find a point in 
this narrow region, one must have T > For clusters 
which are not one-dimensional, so that F' (z) is less sin- 
gular, a smaller value of T may be used. A final check on 
the value of T needed is to look at the number of times 
a trial angle is accepted, as shown in Table I. 

In all simulations we used Aq — 0.5. We chose a — 2/3, 
suggested by Davidovitch et. al. pl| as minimizing nu- 
merical error due to regions in the fjords of the cluster in 
which the derivative of F{z) varies rapidly. An improved 
technique for dealing with these regions involving an ac- 
ceptance window was suggested by Stepanov and Levitov 
|12| . However, for 77 > 1, we automatically avoid regions 
with large |F'(2;)|, so that we do not expect to have to 
worry about this problem; this justifies the very simple 
plotting technique employed to produce the images of 
the clusters below. We chose to plot only the images of 
certain points on the unit circle chosen to lie near sin- 
gularities of the mapping, each point corresponding to a 
single growth step. 

We used a parallel computer to generate many real- 
izations of each cluster. However, generation of a single 
large cluster is a process that parallelizes very efficiently. 
Each node picks random angles and computes the Ja- 
cobian, returning the results to a central process, which 
does the acceptance calculation, and returns the chosen 
angle to the other nodes. This will increase cluster size, 
but at a cost in statistics. 

Results — In Fig. 2, we show a large cluster with rj = 
2, TV = 50000, T = 30 for purposes of illustration. In 
Figs. 3,4 we show clusters with 77 — 3.5,4.5 and N — 
10000, T = 500. A difference in the "fuzziness" of the 
two clusters is clear, with the cluster at 77 = 3.5 having 
many more small side branches. The single side branch 
in Fig. 4 is not the only possibility for this rj. Other 
realizations show either a cluster with no side branches, 
or with a single branching near a growth tip. 

For quantitative results, we have generated 13-14 re- 
alizations of clusters for each 77 = 3,3.5,4,4.5,5, with 
N = 10000, T = 500, and 8 reahzations of clusters for 
?7 = 2 and = 15000, T = 200. provides a measure 
of the linear size of the cluster; it may be shown that the 
radius of the cluster is at most 4 times Fi. We calcu- 
lated the dimension D by averaging Fi over realizations 
of clusters and using linear regression on a log-log plot to 
fit 

Fi("^cxni/^. (4) 



In Fig. 5 we plot Fi and the rms fluctuations in Fi for 
77 = 3 on a logarithmic scale; the rms fluctuations have 
a smaller slope, so the relative fluctuations in Fi tend 
to zero, as found previously for rj — 1 ||ll|. The use of 
-Fi to determine dimension avoids problems with finite 
size effects in box-counting that can lead to a spurious 
dimension less than 1, as found for some t] in previous 
studies iQ. 

While for 77 = 3, Fig. 5 is close to straight, for rj ~ 4 
there is a curvature, as shown in the solid line in Fig. 6 
(the dashed line is a fit discussed below). To quantita- 
tively measure this curvature, indicative of corrections 
to scaling, we performed the linear regression 4 times for 
each 77, first on the full set, then on the last half, quar- 
ter, and eighth of the data set. The results are shown 
in Table II; we also show the rms fluctuations in Fi af- 
ter the 10000-th growth step. For 77 < 4, the corrections 
to scaling are small, and the dimension is relatively con- 
stant, sometimes increasing and sometimes decreasing in 
the last part of the data set; these dimensions are close 
to those found in previously [l^-p^ . 

For 77 > 4, there is a clear trend for the dimension to 
approach unity at longer times, with the trend most clear 
for 7] > 4. We test the trend against the RG prediction for 
the corrections to scaling. The leading irrelevant variable 
in the RG is a tip splitting rate g, which is irrelevant 
for 77 > 4, with scaling dimension (4 — 77)72. At 77 = 4 the 
tip-splitting is marginal, which may be shown to lead to 
Fi = ciTt/ (log (71/02))'^^, where C3 is a universal constant 
determined by the /3-functions and C2 is nonuniversal. 
The dashed line in Fig. 6 shows the flt at 77 = 4 with 
this scaling form; the flt is indistinguishable from the 
data over more than three decades, with only a slight 
difference on the flrst 5 data points out of 10000. (Within 
the RG, the /3-functions are determined numerically, so 
we keep C3 as a fltting parameter; the flt agrees with the 
numerical result from the RG). For rj > A the leading 
correction to scaling is Fi = ci7i(l -I- (n/c2)^'*~''-'^^). A 
flt with this form is shown in Fig. 7 for 77 = 4.5, with 
the flt distinguishable from the data over only roughly 
the flrst 30 data points out of 10000. (We checked the 
RG prediction for the scaling corrections by attempting 
to flt the numerical data with other forms; the results are 
much less accurate). 

Conclusion — We have presented a fast algorithm for 
studying the dielectric breakdown model. Numerical re- 
sults using this algorithm, using the RG corrections to 
scaling, support an upper critical 77c = 4 for which the 
DBM clusters become one-dimensional. The existence 
of flnite size corrections for rj > 4 makes it desirable to 
simulate even larger clusters, but, while we have greatly 
improved the cluster size over previous work, the results 
in this paper are close to the limit of what can presently 
be attained using the conformal mapping techniques. 
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TABLE IL For various values of rj, calculated dimension 
based on different parts of data set; also, rms fluctuations in 
_Fi after last step, normalized by 



V 


full set 


last 1/2 


last 1/4 


last 1/8 


rms 


2 


1.433 


1.4256 


1.435 


1.4522 


.039 


3 


1.263 


1.264 


1.262 


1.243 


.056 


3.5 


1.170 


1.143 


1.155 


1.162 


.078 


4 


1.128 


1.090 


1.078 


1.071 


.072 


4.5 


1.101 


1.090 


1.066 


1.035 


.039 


5 


1.068 


1.030 


1.025 


1.009 


.046 
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TABLE I. For various values of r], N, T, the total number 
of acceptances, and number of growth steps on which at least 
one acceptance occurred, for a typical run. 



V 


N 


T 


Total 


Growth Steps 








Acceptances 


with Acceptance 


2 


15000 


200 


1291738 


15000 


2 


50000 


30 


580805 


49902 


3 


10000 


500 


837291 


10000 


3.5 


10000 


500 


442854 


10000 


4 


10000 


500 


226170 


9990 


4.5 


10000 


500 


143031 


9950 


5 


10000 


500 


99702 


9826 



3 



4000 




FIG. 5. Average of Fi (upper curve) and rms fluctuations 
in Fi (lower curve) as a function of n for = 3. 
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FIG. 2. Gluster with r? = 2, JV = 50000, T = 30. 




FIG. 3. Cluster with r) = 3.5, N = 10000, T = 500. 
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FIG. 7. Plot of Fi versus n for ?7 = 4.5. The dashed hne 
fit including corrections to scaling. 



